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The periodic s-d(f) exchange model is characterized by a wide variety of interesting applications, 
a simple mathematical structure and a limited number of reliable approximations which take care of 
the quantum nature of the participating spins. We suggest the use of a projection-operator method 
for getting information perturbationally, which are not accessible via diagrammatic approaches. In 
this paper we present in particular results beyond perturbation theory, which are obtained such that 
almost all exactly known limiting cases are incorporated correctly. We discuss a variety of possible 
methods and evaluate their consequences for one-particle properties. These considerations serve as 
a guideline for a more effective approach to the model. 



I. INTENTION 

The interplay of localized magnetic moments (of d- or 
f-type) with itinerant s-electrons of a partially filled con- 
duction band is of indisputable importance for the expla- 
nation of many effects in condensed matter physicsiSi. 
In recent years and in the context of diluted magnetic 
semiconductors it even became a driving mechanism for 
electronic applications (spintronics)^. For a theoretical 
description the Hamiltonian 

H^no + nsf^Y.J2 ^Is- -jY^ct^-s.. (1) 

is probably the simplest choice possible. The local- 
ized magnetic moments are represented by quantum- 
mechanical spin operators and interact with the spin 
of the conduction electrons cr.j via an exchange interac- 
tion J. The dispersion of the latter is determined by 
the hopping integrals T^- . In recent years one has referred 
to the Hamiltonian as the Kondo-lattice model, nev- 
ertheless sf-, sd- or (for very large J) double-exchange 
model are probably more suitable descriptions. 

Despite the large variety of systems to which the s- 
f model can be applied^, a convincing approximation 
scheme is still missing. Previous attempts developed 
along two principle directions: On the one hand the 
Green's function hierarchy of equations of motions has 
been restricted by some decoupling schemes^. Even 
though these efforts allowed a fair description of mag- 
netic semiconductors such as europium chalcogenidegi, 
this kind of approximation always suffers from a lack of 
controllability. On the other hand, a classical treatment 
of the localized spins lead to some substantial results 
obtained using dynamical mean field theory^ or Monte 
Carlo techniques^. However, it has been shown that the 
quantum-mechanical character of the spins has indeed 
a substantial impact on the electronic properties of the 
concerned materialsiSiii. 

Therefore, an analytical and numerical treatment of 
the s-f model should be aspired to, which retains the 
quantum nature of the spins and ensures the control- 
lability. A perturbational approach would be a good 
candidate for the latter. However, the incorporation of 



quantum-spin operators causes difficulties since Wick's 
theorem, which is generally used to evaluate Feynman 
diagrams, is not applicable to this class of operators. 
This is probably the reason why up to now only few at- 
tempts exist in this direction. These are mostly limited 
to special situations like one-dimensional systems or half- 
fillingi2.i3. 

We suggest to circumvent these difficulties by using 
the projection-operator method (POM) as introduced by 
MoriMii^. Its goal is the expansion of resolvents such 
as the one-particle Green's functions in a continued frac- 
tion. The further one goes in the expansion of the con- 
tinued fraction, the higher the accuracy of the results 
obtained. It has already been applied successfully for 
a weak-coupling approach to the Hubbard model where 
it leads to convincing resultsiSiii. Its application to 
the s-f model would be an obvious development, even 
though straight forward test calculations were not very 
successful 

A weak-coupling approach to the s-f model is undoubt- 
edly very interesting, both from an experimental and 
a theoretical point of view. First of all, the parameter 
regime is probably important for diluted magnetic semi- 
conductors. Secondly, the well-known RKKY exchange 
mechanism.^®'^^'^°, usually used for a qualitative descrip- 
tion of magnetism in these materials, is also nothing else 
but the application of second order perturbation theory 
to Hamiltonian However, the goal of this approach is 
the determination of the ground state energy, leading to 
an effective Heisenberg interaction. In contrast to that, 
here we will show that we are able to analyze the much 
richer spectrum of dynamical properties of the s-f model 
to the same order of the coupling constant and beyond 
it. 

Whenever a second order perturbation theory (SOPT) 
description is performed in many-body theory there are 
principally three different ways of treating emerging 
propagators: a conventional SOPT uses only free prop- 
agators, a SOPT relative to Hartree-Fock (HF) replaces 
these propagators by the corresponding mean-field ex- 
pressions and a self- consistent treatment of the SOPT 
replaces all propagators by (functionals of) the full prop- 
agators as obtained in the previous step of iteration. A 
priori it is usually not known which version yields the 
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most reliable results. Of course, the self-consistent ver- 
sion includes a summation of more diagrams than the 
other methods. However, since only a partial class of 
diagrams is summed, it is unclear which important dia- 
grams are being missed out or canceled, double-counted 
or even taken wrongly. There are more profound consid- 
erations as put forward by Kadanoff and Baym^^ argu- 
ing that the self-consistent approach is a conserving ap- 
proximation which automatically satisfies the Luttinger 
theoremjS? and Fermi-liquid properties. However, it is 
uncertain to which degree this applies to a model which 
is not exclusively composed of fermions. 

Even for more established models such as the peri- 
odic Anderson model (PAM), the Falicov-Kimball model 
(FKM) or the Hubbard model (HM) the question of 
the most appropriate version has been discussed inten- 
sively. For the PAM Yamada and YosidaiSSiSlS^ started 
the perturbational investigations directly by considering 
deviations from the non-magnetic HF solution. Later 
Schweitzer and Czycholl^^ were able to numerically com- 
pare this approach with a self-consistent SOPT. Even 
though this version obeys more of the Luttinger sum 
TulesS^, the self-consistent version (in contrast to the ver- 
sion relative to HF) failed to show the one-particle peaks 
near Ef and Ef + U in the /-electron spectral function. 
The FKM is another example of a model where a not 
fully self-consistent SOPT treatment qualitatively repro- 
duces exact results, whereas the self-consistent SOPT 
does not,^^ For the HM the situation is slightly more 
complicated. One can show that SOPT relative to HF 
does not yield a metal-insulator transition and does not 
show a breakdown of the Fermi-liquid behavior. On 
the other hand, a straightforward application of a self- 
consistent SOPT does not reproduce the Hubbard bands 
in the atomic limit More sophisticated methods such 
as the interpolation scheme of Edwards and Herta2L2S 
(a version relative to HF) or the iterative perturbation 
theory (IPT) of Georges and Kotliar— (a self-consistent 
version) are required. For the latter approach the re- 
striction to half-filling has been removed by the modified 
perturbation theory (MPT) of Kajueter, Kotliar'^^ and 
Potthoff, Wegner, NoltingSi. The MPT is probably the 
most convincing analytical approach to the HM.^- 

With the present paper we extend this kind of discus- 
sion to the electronic part of the periodic s-d(f) exchange 
model. We will argue that it is indeed a self-consistent 
ansatz for the electronic self-energy which is the most 
promising for this model. Other possible weak-coupling 
approaches are ruled out after a direct comparison with 
the results of our method of choice. It can be shown 
that minor changes in the analytical method have dras- 
tic effects on one-particle properties, such as the density 
of states. We believe that a more profound analysis of 
the s-f model (e.g. by a combination with band-structure 
calculations) can be based on these considerations. 

The article is organized as follows: As a starting point 
we derive in section ^ an SOPT for the s-f model which 
makes use of the POM. In section IIIII we will study in 



some detail the exactly soluble limit of a single conduc- 
tion electron in a ferromagnetically saturated semicon- 
ductor. This limit is an excellent testing ground for the 
implementation of the POM. Even more important, the 
experience with other models shows that it is indispens- 
able to have non-perturbative, exact statements which 
can be used to judge the quality of the results obtained. 
In a next step (section IIV|) the experiences for an im- 
provement of this limit are generalized to arbitrary pa- 
rameter configurations. With all these preparations we 
are able to present the results of the self-consistent calcu- 
lation and a comparison with other methods in sections 
IVl and IVII The findings are summarized in section IVIII 

It is worth mentioning, that our approach seems to be 
related to a moment-conserving interpolation scheme of 
the self-energy as published by Nolting et al.^'^'^. There, 
a general structure of the electronic self-energy, which 
looks similar to the one presented here, has been found 
by systematically studying all known exact statements 
on the s-f model. However, their analysis is focused on 
the low-density limit and ensures the correctness of these 
statements for n — > (or n — > 2) only. In contrast to 
that the approach given in this paper concentrates on 
the weak-coupling behavior. Indeed, independent of the 
occupation number the correctness of the self-energy up 
to order in the coupling parameter and up to order 
E~^ in the high-energy expansion is guaranteed! Addi- 
tionally, we can fulfill the same criteria for n — > as 
given in the above mentioned publication. Nevertheless, 
the two approaches are not identical even for n = 0, but 
otherwise arbitrary parameters. 



II. SECOND ORDER PERTURBATION 
THEORY 

As mentioned above, we use the projection-operator 
method (POM)^'*'^^ since it allows an expansion of resol- 
vents without the use of Wick's theorem. The approxi- 
mation consists in considering only a physically relevant 
subspace of the Liouville space. With the simplest choice 
the Liouville space is spanned by single-particle states 
cj^^). Accordingly the projection operator and its or- 
thogonal complement are defined as 



P 



L) Ua and Q^l - cL) (cL 



(2) 



These definitions require the existence of a scalar prod- 
uct, which in our calculations is conveniently chosen to 
be the thermodynamic average: 



iA\B)^([A+;B] 



(3) 



Within the POM the one-particle Green's function is 
given by the following dynamical equation 
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[r!k. + Mk.(c^)]XkJ' 
(4) 
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where uj = E + iO"*" and the Liouville operator C with its 
property C\A) = \[H,A]_) has been incorporated. For 
the choice ^ the susceptibihty matrix xwa = {(^ka\^ka) 
is particularly simple: Xkcr = 1- The frequency matrix 



C 



(5) 



on the other hand, corresponds to the mean-field result 

1 



(6) 



for the Green's function, 
ations J = hJ and 



Here, we have used the abbrevi- 
2^ c^.i^ = ±1. All the interesting physics 
is included in the memory matrix 



M, 
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oj — QCQ 



(7) 



which again has the structure of a resolvent, resulting in 
a form of Gko- in Q involving continued fractions. 

The expression for the memory matrix cannot be 
treated exactly. However, at this stage we are only aim- 
ing at a perturbational expansion of the self-energy Sko- 
in the form 



(8) 



This allows some simplifications. In (UJ already the Hsf- 
contribution in IQ-Ccj^^) gives rise to a factor J^. Hence, 
any approximation of the Liouville operator in the de- 
nominator is still correct in this order and is thus consis- 
tent with (jSJ. A conventional SOFT implies a replace- 
ment of C by its free part Cq. A SOFT relative to HF 
is given by a Liouville operator that corresponds to the 
Hamiltonian 
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It turns out that this procedure is equivalent to a replace- 
ment of the Green's functions at the right side of (|1U|) by 
the full propagators as obtained in the previous iteration 
cycle. In we use the additional approximation that 
the memory-matrix is summed over k and hence only a 
local self-energy is considered. 



III. AN EXACT SOLUTION 

The model can be restricted to the limit of a ferromag- 
netically saturated semiconductor. This limit is char- 
acterized by two mathematical consequences: Firstly, a 
semiconductor is defined by an empty conduction band 



at zero temperature, hence 



) = 0. Secondly, fer- 



romagnetic saturation leads to trivial spin expectation 
values: (. .. S^) = NS {...),{.. . S+) = 0. An applica- 
tion of these simplifications to the result (|10|l yields the 
self-energy 



Sko- — So- — 



(12) 



One can see directly that in this limit a self- consistency 
iteration does not yield any further results. This is 
because the a =| Green's function is uniquely deter- 
mined by a =t propagators, which have the iteration- 
independent self-energy Sj = ~JS. 

Due to the restrictions of this limit the memory ma- 
trix can in fact be treated more accurately. Following 
the intention of continued-fraction expansion the mem- 
ory matrix can itself be considered as a resolvent, to 
which the concept of the POM is applied. 
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In both cases we obtain a similar result 

1 



7k. = -(5^)^4r"'+^E(^v^><C(io) 



(0/MF) 
k+q,-o-' 



where = 



izcr5^. For the expectation values con- 
tained in (|10|l we make use of the fact that we aim for 
a result correct to second order in J and evaluate them 
using the eigenstates of the free/mean- field system. 

A self-consistent SOFT on the other hand can be ob- 
tained in the same manner as suggested by Bulk and 
Jelitto^ for the Hubbard model. Within this procedure 
the unperturbed part is altered in each iteration cycle by 
the memory-matrix of the previous cycle: 



(N+l) 




(N) 



kcr- 



(11) 



The higher-order memory matrix M will have a form sim- 
ilar to that given in ijJJ). Again the Liouville operator C 
in the denominator should be approximated to allow for 
an analytical solution of the associated geometric series. 
According to conventional perturbation theory it is re- 
placed by Cq, the action of a free, undisturbed system of 
electrons. After sophisticated calculations, which will be 
published elsewhere, the self-energy is obtained as 



i--/(i-^)^E<^(' 



(14) 



This is certainly an improvement of 1)12(1 and contains the 
previous result if expanded up to order J^. 

Even though we called it perturbation theory, it is how- 
ever not correct for the next order in J. The exact J^- 
contribution to the self-energy can actually be shown to 
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be 



S ^ — T 

q 



[g(°)(c.) 



(15) 

As a matter of fact, the second sum in (|15|l is a diverging 
contribution. This ah'eady becomes apparent if one looks 
at its imaginary part, rewrites the k-sum as an integral 
over the free DOS, separates a Lorentzian and considers 
the fact that 0+ is infinitesimally small. We were able 
to show that after a summation over all orders in J the 
diverging terms cancel. Nevertheless, (|15|) demonstrates 
that for the s-f model a strict perturbation theory is only 
possible up to second order in J. The different orders of 
an expansion of the continued fraction within the POM 
apparently do not have this limitation. 

The result (O can be further improved, if a perturba- 
tion theory relative to HF is chosen. Without going into 
the details we just provide the result as: 



-J ZaS + 5al 



(16) 



This expression is actually identical to the result of an 
exact calculation, where the Liouville operator has not 
been reduced or altered. Within the POM we were able 
to perform the derivation of the memory matrix ((T)) us- 
ing the complete operator C = Co + Csi- However, it 
is not necessary to give the lengthy calculations here 
since its result ()16f) has already been verified by other 

methodS|35,36i37,38,39 

The self-energy (|16f) corresponds to an exact eigenstate 
of the Hamiltonian . For the spin-down electrons this 
eigenstate, which is the ground stated for antiferromag- 
netic couphng (J < 0), is called magnetic polaron^i. Its 
interesting and nontrivial dynamical features, which give 
rise to a scattering part and a quasi-particle peak in the 
density of states have been discussed in detail by Nolt- 
ing et ali^. Apparently, it is possible to retrieve these 
features within the projection-operator formalism. The 
reason why already an approximation yields the correct 
result is the fact that the resulting two-dimensional Li- 
ouville subspace is sufficient to completely describe the 
physics of a ferromagnetically saturated semiconductor. 

Now we have obtained several approximate forms of 
the electronic self-energy. Formula H12|l provides an ex- 
pression for the first and second order in the coupling con- 
stant J . The (diverging) third order is given in H15|) . An 
improvement of the SOPT is given in H14|l on a conven- 
tional way and in relative to HF. One can compare 
these self-energies by looking at their quasiparticle den- 
sities of states (QDOS). Obviously, the spin-t spectrum 
is always a mean-field-shifted free DOS. Hence, we can 
limit ourselves to the spin- J, spectrum, which is shown in 
figure ni 

In this figure a relatively large J — 2 J has been chosen 
to reveal the differences more clearly. If one compares the 
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FIG. 1; Comparison of the J,-QDOS for different self-energies 
in the limit n = and (5^) — S. The dotted and dashed 
line are the results of SOPT (I12L resp. conventional and 
relative to HF. The dashed-dotted and the solid line give the 
corresponding results H14^ / l|16|l for the next step of the POM. 



DOS of the conventional SOPT and the SOPT relative 
to HF with that of the exact solution, one gets the im- 
pression that the first one is the better approximation. 
However, if one compares these two approaches for the 
second step of the POM (equivalent to a larger relevant 
Liouville subspace) it is clear that the version relative to 
HF ha,s to be preferred, since only this one gives the exact 
result. Nevertheless, result (|14() is already a satisfactory 
approximation. As mentioned above a self- consistent cal- 
culation is redundant for the discussed limit. 



IV. GENERALIZATION TO ARBITRARY 
CONFIGURATIONS 

After these considerations on the magnetic polaron we 
return to the discussion of the second order perturbation 
theory (SOPT) given in formula (|10|l . As argued before, 
the limit of a ferromagnetically saturated semiconductor 
can be used to check the quality of this result. In this 
limit the SOPT-result (|12|l turns out to be only a poor 
approximation of the exact solution as demonstrated in 
figure n Since the magnetic polaron is indeed an impor- 
tant feature of the s-f model, ways of improving (|8I10() 
should be considered. 

In the previous section we explained how an improve- 
ment can be achieved within this particular limit. A 
proper application of an additional step within the pro- 
jection operator method finally leads to the expressions 
(|14|l and We generalize the analytical structure of 

these results to the following ansatz for arbitrary band 
occupations: 



1 - bkall^aiE)' 



(17) 
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Although obtained in a completely different manner, this 
is exactly the kind of a modified perturbation theorySi 
(MPT) which has turned out to be the most promising 
analytical approach to the Hubbard model'^^ . 

For the s-f model it has two advantages: First, it does 
not destroy the correctness of the second order term pro- 
portional to J^, but gives the freedom to fit the pa- 
rameters ak(T and ^ko- such that further criteria are ful- 
filled. Secondly, since the SOPT-result for the self-energy 
© automatically reproduces the first three moments of 
the corresponding Green's function correctly, the choice 
flko- = 1 will ensure the same for the ansatz (|17|l . 

It remains to determine the parameter foko-- The most 
straightforward choice merely ensures the correctness of 
the ferromagnetically saturated semiconductor ( &kcr = 
J/2S ). However, the resulting densities of states have 
unphysical features. Furthermore, we have learned from 
the Hubbard model=^ that a fit to the spectral moments 
of the Green's function yields more promising results. A 
6k(T which is determined by 



[qCcI^\qCcI^ 



(18) 

ensures that the third coefficient of the high-energy ex- 
pansion of H18|l is identical to the one of the exact self- 
energy. As explained in more detail in appendix 1X1 this 
fit is correct for the first four moments of the Green's 
function. 

In an MPT which is based on a conventional SOPT the 
Liouville operator Cq in H18|l is understood to correspond 
to the free part of the Hamiltonian Hq. Then only the 
J'^-contribution in the numerator of feko- remains to be 
evaluated and one obtains 



SiS + l)~z„ (5^) - {S'-y {z„ {S') + l) + q^ 



SiS+l)-z, {S^)-{S^Y + 2p. 



(19) 

where pa and qa are sets of further correlation functions, 
which are given in appendix ^ but have the property 
to vanish in the limit n — > 0. It is instructive to study 
this limit in more detail. On the one hand it can be 
combined with the additional constraint of ferromagnetic 
saturation (S"^) = S. If the obtained 6| is placed into the 
MPT- ansate (|17() . then the self-energy becomes identical 
to the one given in 1)14(1 . On the other hand, one can 
consider the zero-bandwidth situation ek = Tq within 
the limit n ^ 0. For a dispersion-less Green's function 
all summations in expression (|10|l for 7kCT can readily be 
performed and in this limit the self-energy becomes 

,2 [^(^ + l)-^.(^^)-(^^)']^ 



1- J(z^(5-) + l)^ 



(20) 

Comparing this form of Sg. to the result of an exact cal- 
culation available for the zero-bandwidth limit^'-- reveals 



that the expression is correct. 

For a SOPT relative to HF the expression for Cq in 
equation ((18|l contains an additional term according to 
the mean- field contribution Eq- = —Jz^ (S^) in J^. This 
yields a correction J&o- — 

[ ((5^)') - {S'f ] + (S-^S^) + 2z„ (S^) 



S{S+l)~z^{S^)-{S-f + 2p, 

(21) 

to the former result (|19() . We will focus again on the limit 
n — > 0. If combined with the additional constraint of fer- 
romagnetic saturation, the exact result ()16|l is obtained. 
A consideration of the zero-bandwidth limit yields the 
following expression for the self-energy 

.Py 

\2i 1 



(22) 



[S{S+l)-z,{S^')-{S^r 
J2 



1 - J 



X 



Apart from the correction terms X and Y this is again the 
exact result (120(1 . However, both terms are proportional 
to expressions, which vanish in the paramagnetic regime. 
Therefore, they also vanish for the zero-bandwidth limit 
for which the assumption of any finite magnetization does 
not lead to consistent results^i 

In the discussion so far we have tested our MPT- ansatz 
((17(1 in the limit n 0. One can repeat the same trans- 
formations for the opposite case n — > 2. By doing this 
one will notice, that the same formulae are obtained. 
The only difference is the change of the sign of a and 
of 6k(T- This is due to particle-hole symmetry in the sys- 
tem. Therefore, in the same sense as for n = our MPT- 
ansatz ((17(1 fulfills the limit of the magnetic polaron and 
the zero-bandwidth limit for n = 2. 



V. SELF-CONSISTENT RESULTS 

Using the MPT approach of the previous section we 
have the possibility to generalize the improvement to 
the POM for the limit of a ferromagnetically saturated 
semiconductor to arbitrary parameter regimes. In sec- 
tion mil we argued that a self- consistent calculation is 
redundant in this limit. This does not hold for the gen- 
eralized version. Here 7ko-5 as given in ((10(1 . does not 
vanish for spin— "f electrons and consists of propagators 
for both kinds of spin-directions. Consequently, we have 
performed a self- consistent immeTical iteration of the self- 
energy. This was carried out along the lines sketched at 
the end of section ^ Additionally, the MPT-parameter 
ba- has to be adjusted such that So- in ((21(1 describes the 
full self-energy and not only its mean-field part. 

The details of this procedure are discussed in section 
IVII If properly performed it is an "upgrading" of the 
perturbation theory relative to HF in the sense that its 
properties are maintained. In particular, the evalua- 
tion for the ferromagnetically saturated semiconductor 
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FIG. 2: Dependence on band occupation n for a self- 
consistent MPT. For free electrons a sc-DOS with bandwidth 
W is chosen. The parameters are as given. The vertical lines 
indicate the positions of the Fermi energy, to allow for the 
different values of the band occupation. 
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FIG. 3: Change of the QDOS with the magnetization (S^) 
for a self- consistent MPT. The position of the Fermi energy 
for the highest and lowest (5*^) is marked by vertical lines. 



yields the previous and exact results of the magnetic po- 
laron. Additionally, the atomic limit is fulfilled for the 
empty and completely filled conduction band, and the 
particle-hole symmetry of the system is conserved by the 
ansatz. In other regimes it is ensured that all results 
are correct at least to order J^. However, because of 
the self-consistency the method incorporates more cor- 
relations and scattering effects than a straightforward 
second-order perturbation theory description does. For 
all these reasons we believe that the self- consistent MPT 
is not only correct in the weak-coupling regime, but also 
for moderate values of J. 

We discuss the results in terms of the quasiparticle 
densities of states (QDOS). According to our ansatz the 
QDOS is correct for n = 0, {S^) = S, however its vari- 
ation with a change of these parameters is of particu- 
lar interest. Figure |21 shows the dependence on the first 
parameter. The dependence on the magnetization (con- 
nected to temperature via a Brillouin function) is given 
in figure 121 In both cases a medium value has been cho- 
sen for the fixed parameter, respectively. It goes without 
saying that our calculations are also thermodynamically 
self-consistent. The iteration ensures that the values of 
the correlation functions are consistent with the obtained 
one-particle Green's function. Additionally, the chemical 
potential is adjusted to the desired particle number. Its 
position is indicated by vertical lines in the figures. 

In particular the change of n in figure El has remark- 
able consequences for the QDOS. The clear dependence 
on the filling of the conduction band points out strong 
correlation effects, induced by the coupling J. For n = 
the structure of the QDOS is closely related to the ferro- 
magnetically saturated semiconductor. In particular the 
spin-t spectrum has the shape of the free (simple cubic) 
DOS and the scattering part of the spin-| spectrum can 
clearly be seen. Only the polaron subband shows a de- 
formation, due to finite-lifetime effects. Excited spin-| 



electrons can enter the energy region of the polaron, flip 
their spin and absorb a magnon since we are not close to 
saturation. 

If the chemical potential (and accordingly the band 
occupation) is increased, the spectral weight is redis- 
tributed between both subbands. For the chosen set 
of parameters the changes with n are most noticable in 
the spin-| QDOS, where the upper subband steadily in- 
creases in importance at the expense of the lower sub- 
band. A sharp jump in the QDOS close to the pseudo- 
gap remains a striking feature for all values of n. It is also 
interesting to note, that the lower band edge is shifted 
by some 0.1 eV in the spin-t QDOS, whereas it remains 
at almost the same position for the spin- J, QDOS. This 
behavior is very much different in an MPT relative to 
HF^^ and is a hint that in the self- consistent MPT it are 
mainly the majority-spin electrons that experience strong 
correlations. 

As the band occupation approaches half filling {n = 1), 
the point-symmetric form of the QDOS nicely represents 
the particle-hole symmetry of the system. The character 
of the upper spin-| subband becomes identical to the 
lower spin-f subband, since the latter is the polaron band 
for n — 2. For the same reason we skip the plots for 
ri > 1, they can be obtained from the band occupations 
2 - n. 

A higher value for the coupling strength J is chosen 
in figure (31 Therefore, the scattering and the polaron 
subband are well separated already for a nearly satu- 
rated system. The gap remains present for all tempera- 
tures employed in the calculations. There are only small 
changes of the position of the bands as a function of the 
magnetization. Nevertheless, the edge of the lower spin- 
t subband shifts to lower energies if the temperatures is 
lowered from T = Tc ((5^) = 0) to smaller values T ^ 
(maximum (S^)). For semiconductors such an effect is 
known as the red shift of the optical absorption edge. 
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In metals, since the lower spin-| subband is shifted in 
the opposite direction, it leads to a polarization of the 
conduction electrons of over 60%. The existence of en- 
ergy regions well below the Fermi edge occupied entirely 
by majority-spin electrons is a remarkable result. Similar 
effects have also been reported in other approximations"''^ 
when studying half-metallic ferromagnets. However, the 
continuous shift of the chemical potential with magneti- 
zation prevents 100% polarization of the conduction elec- 
trons in our calculations. The dependence on the chemi- 
cal potential is such that the effect disappears completely 
for smaller values {n 0), where the t-QDOS and the 
j-QDOS occupy the same energy region (see figure ISJ. 

The situation at the lower edge of the upper subband 
in figure 13 is less systematic. An extra peak in the spin- 
t QDOS obtains its maximum for (S*^) « 0.3 and then 
vanishes again. The physical interpretation of this fea- 
ture is not yet clear. In this context we also have to 
mention that the choice of the parameter with its 
explicit dependence on the self-energy in (|21|l substan- 
tially increases the numerical effort. In the energy region 
just mentioned it is particularly difficult to obtain con- 
vergence. 

The distribution of the spectral weight, on the other 
hand, is indisputable. It shows the transition from the 
ferromagnetically saturated configuration (dashed lines), 
which even for n — 0.4 clearly displays the features of the 
exact solution in section ITTTl to the paramagnetic regime 
(solid lines), which has to be symmetric with respect to 
the X-axis. Again the more profound changes are ob- 
served for the majority-spin electrons. The increasing 
spectral weight of the upper spin-f subband can be ex- 
plained with higher magnon-numbers in this regime. 

An artifact of our method is the fact that for both spin 
directions two subbands are always obtained. With other 
approaches24, one sometimes observes a third band. This 
is explained by atomic- limit calculations, where for finite 
band-occupations always three out of four subbands have 
non- vanishing spectral weight. It needs further modifi- 
cations of our method to retain these features. At the 
present stage the atomic limit is only correct for n = 
and n — 2. 



VI. THE PROPER METHOD 

The last point brings us to an assessment of our self- 
consistent approach. Apart from the above mentioned 
catalogue of analytical properties a comparison with 
other conceivable weak-coupling approaches is desirable. 

A comparison with an MPT which uses 7kcr obtained 
by conventional SOPT is straightforward. From the ana- 
lytical considerations in section Hvl we can conclude that 
even an MPT based on an SOPT relative to HF should 
be preferred as compared to one based on a conventional 
SOPT. This is because the former correctly incorporates 
the important limiting case of the ferromagnetically sat- 
urated semiconductor, whereas the latter does only re- 



produce, in this limit, the less accurate expression H14|l . 
The discussion of the atomic limit does not provide an 
argument in favour of one the approaches, since at the 
end both yield physical expressions of the same quality. 
In these limits our self-consistent approach has the same 
properties as the MPT relative to HF. 





<S"> = 0.1) 




<S'> = 0.1 




<S"> = 0.2 




<S') = 0.3 




- - - (S'> = 0.4 




<S'> = 0.5 




FIG. 4: Change of the QDOS with the magnetization (S^). 
The MPT-calculations are based on a SOPT relative to HF. 
The parameters are chosen as in figure |3 

Our comparison with an MPT relative to HF is based 
on numerical results with this method. For the same pa- 
rameters as in figure 131 we obtain the set of QDOS given 
in figure ^ Even though its main features look sound 
again, there are a set of minor aspects which make this 
approach questionable. The most obvious drawback are 
prominent and unexplainable peaks close to the Fermi 
energy. The functional dependence is not smooth and 
dominated by the free DOS. Also for magnetizations be- 
low saturation a gap in the spin-down QDOS is expected 
from other theoricsMi, but this does not exist in this 
approach. Additionally it is noteworthy that the onset 
of the spin-up QDOS starts for smaller energies as com- 
pared to that for the spin-down QDOS. 

In the self- consistent MPT of the previous section most 
of these peculiarities are not present. A comparison with 
figure 131 permits the conclusion that the iteration of the 
self-energy yields shapes of the QDOS which are broader 
and smoother. This can be understood analytically and 
is closely related to the fact that the self-energy be- 
comes complex by iteration. These findings make the 
self- consistent version more reliable. 

However, there are several possibilities to incorporate 
self-consistency into the MPT. An ambiguous point is the 
order of the applied steps. In contrast to the calculations 
of the previous section one could start with the expres- 
sion of the MPT as obtained in an approach relative to 
HF and continue by dressing all included propagators as 
full Green's functions. Hence, the parameters b^^ are 
fixed, which is still correct for a fit to the high-energy 
expansion up to order E~'^. The consequences for the 
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J= 1.0 eV 
;/' W=1.0eV 

S = '/2 

n = 0.4 



FIG. 5: Dependence on magnetization (S^) for an MPT based 
on a self-consistent SOPT. The parameters are as in figure|21 



densities of states are shown in figure [S] They look rea- 
sonable for values of (S*^) close to saturation. However, 
the single, broad, elliptic band which emerges close the 
paramagnetic regime is a surprising feature. Not only 
is this result inconsistent with other approximations of 
the s-f mode&Si2i24, moreover it is incorrect in the zero- 
bandwidth limit. This is because the same structure re- 
mains in the limit W —>■ 0/n —>■ 0, where two narrow 
bands around —JS and J{S+l) are expected as indicated 
by exact calculations''*'. The shortcoming can already be 
seen analytically when looking at equations (|20|l or (|22|l . 
If the propagators are dressed without a change of ba-, 
then an equality is not possible with the exact solution: 



1 



' E-Tn-S , T9 ' 



1 - J 
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B-Tn 



1 - J 



(23) 



This conflict can be resolved if the order of arguments 
is changed. Now the starting point is a self-consistent 
SOPT and only afterwards the result is fitted to the high- 
energy expansion. The consequence is not a higher accu- 
racy in powers of but an additionally dressed fitting 
parameter hu. According to the correction given in l|21() 
it now contains the full self-energy and not only its mean- 
field contribution. This leads to the correct result for the 
zero-bandwidth limit at n = 0: 



l-(J-S) 



1) 



E^Tn 



E-To-T: 



1 - J 



(24) 



E-To 



Additionally, it has serious consequences for all other 
parameter regions. We compare numerical results for the 
three different approaches mentioned above in figure 
The set of parameters W, J, 5" is chosen to be close to 
the atomic limit. The choice n = 0, (S^) = ensures 
exactness in this limit. If the input 7kCT to the MPT 
((T7|l is the result of the SOPT relative to HF (dashed 



n = 0.0 



O 2 
Q 




SOPT-HF 

sc SOPT 

sc MPT 



n 


= 0.0 


n 


= 0.2 


n 


= 0.4 




FIG. 6: Comparison of different MPT approaches for parame- 
ters close to the atomic limit (bandwidth W = 0.5 eV, J — 0.6 
eV, 5* = 3/2, (S"^) = 0.0). (a) dashed line: the MPT is based 
on a SOPT relative to HF, dotted line: only the SOPT input 
is treated sdf-consistently and the parameter ba is not altered, 
solid line: the self-consistency has also consequences for ba- 
For the latter there are two versions to treat the correlation 
functions (CF) in ba- either self-consistent (b) or mean-field 
like (c). The vertical lines indicate the energy-positions of 
the maximums of the subbands. 



line) then two nearly free subbands at the correct posi- 
tions are obtained. There seems to be a third flat band 
between them. To dress only the propagators and not 
the fitting parameter fog- (dotted line) is definitely wrong. 
However, one can convince oneself, that an inclusion of 
the full self-energy in the calculation of (solid line) 
really yields a considerable improvement for n = 0. Here 
one again observes the two narrow subbands, the excita- 
tion energies for an electron that aligns its spin parallel 
{—JS) or antiparallel {J{S -\- 1)) to the localized spin. 

Nevertheless, there remains an uncertainty in the de- 
termination of 6o- as far as the correlation function p„ 
in equation ((T^ and is concerned. Its definition is 
chosen such that it can be calculated with the help of the 
one-particle Green's function Gk^'- 



1 



P-a 



irhJN V 

k 



oo 



dE. 



(25) 



This form of the spectral theorem has the handicap to 
be only applicable for determining a sum of correlation 
functions (|B1|) . Whenever being confronted with a single 
correlation function contribuation to (|25|) . we are forced 
to perform a mean-field decoupling: (S^n^^) = (S^) ■ 
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(^n_^y This is still compatible with the expansion in J 
and E~^. It is therefore straightforward to treat those 
correlation functions for which relations such as H25|) exist 
as accurate as possible, and perform approximations for 
the remaining correlation functions. 

The consequences of this methodology for the atomic 
limit are shown in figure Eb- The QDOS looks sound 
for n = 0.0, but shows a broad non-quasiparticle struc- 
ture between the subbands for n > 0. Its spectral weight 
increases with band occupation at the latter's expense. 
This dependence is qualitatively different compared to 
that of a third, intermediate band in the QDOS for the 
SOPT relative to HF, which remains small for all values 
of n. Here, already for a band occupation of n = 0.4 the 
gap is completely filled. It is difficult to find a physical 
explanation for such a behavior. Small satellite peaks 
between the subbands were also reported for other ap- 
proximation methods^ and were attributed to un-trapped 
electrons which experience the global magnetization (S^) 
as an effective quantization axis. A shift of the spectral 
weight within the intermediate structure as a function of 
the net magnetization has also been observed in our cal- 
culations. However, the missing symmetry in the para- 
magnetic regime and the strong dependence on the band 
occupation does not fit into this picture. The accompa- 
nying shift of the two subbands is also surprising. For 
these reasons we believe that these features are an ar- 
tifact of the approximations used as we know that the 
atomic limit is only correctly included for n = 0. 

In contrast to the methodology to determine be as ac- 
curately as possible, it has apparently a much higher 
priority to treat all included correlation functions on an 
equal footing. Figure shows results for the QDOS 
with the same set of parameters as in figure |BJd. The 
only modification in the theory is a mean-field decou- 
pling of all correlation functions in ba- The effects on the 
QDOS is dramatic, as the intermediate structure com- 
pletely vanishes now. Additionally, the shift of the two 
subbands happens in a comprehensible way: Due to the 
particle-hole symmetry we expect for n — 2 two subbands 
at positions —J{S + 1) and +JS. Since our approach 
apparently only allows for a single bandgap, the change 
from an n = to an n = 2 configuration can only be 
implemented by the system if the two peaks move conti- 
nously to their new positions. Accordingly their positions 
at half filling (n = 1) have to be E± = ±J{2S + l)/2, as 
seen in the figure. Also the redistribution of the spectral 
weight takes place along these lines. 

Based on our experience with the MPT we draw the 
following conclusion: The most promising weak-coupling 
approach to the periodic s-d(f) exchange model is a self- 
consistent MPT. Expression H1U|I dressed with full propa- 
gators GicCT should be used as the input from second-order 
perturbation theory. Only afterwards the parameters Oko- 
and 6ko- in the MPT-ansatz H17|l should be determined 
such that the high-energy expansion is fulfilled to power 
E~'^, which implies that bg- carries a dependence on the 
full self-energy. For the correlation functions entering 



6cr it is important that they are all treated on the same 
footing. This is, at this stage, only possible by using a 
mean-field decoupling. 



VII. SUMMARY AND OUTLOOK 

Within the presented work we have demonstrated how 
the projection operator method can be exploited to find 
an analytical approach to the periodic s-d(f) exchange 
model and that it is indeed a valuable tool in this con- 
text. Nevertheless, we argued that the second order per- 
turbation is insufficient and suggested an improvement 
in the form of an MPT. The principle structure of this 
ansatz results from a study of the limit of the ferromag- 
netically saturated semiconductor. We showed that the 
calculations have to be performed self- consistently. In 
a further step a clarification of the proper treatment of 
the fitting parameter b^ was necessary. At the end we 
were able to make an informed statement, which of all 
possible approaches is the most reliable one. On the one 
hand it is satisfactory that a certain approach was able to 
produce considerably better results than other attempts. 
On the other hand, the high sensitivity of the QDOS to 
the methodology used to treat the correlation functions 
includes the danger of arbitrariness. 

All of the suggested approaches have in common that 
they are correct up to second order in the coupling pa- 
rameter J. Nevertheless, our improvements are non- 
perturbative in the sense that no Taylor expansion is pro- 
vided. On the one hand this is for technical reasons, since 
the absence of a Wick's theorem for spin operators signif- 
icantly complicates the calculation of Feynman diagrams. 
On the other hand already Shastry and Mattis^^ argued 
that a perturbation theory in J would fail because of the 
discontinuities in the physical properties of the model as 
J changes sign. Additionally, we have pointed out that 
the J^-contribution to the exact self-energy in the limit 
of a ferromagnetically saturated semiconductor diverges 
although the sum over all orders yields a finite result. 

The qualitative properties of the densities of states 
presented here are very similar to the findings of other 
approaches^. Since the former results were based on de- 
coupling schemes for Green's functions the approxima- 
tions incorporated into these calculations are difficult to 
control in their quality. With our results we can confirm 
a posteriori and justify these findings. This includes a 
complete set of strong correlation effects discussed there. 

However, in its present state the documented method 
is only an approach to the electronic part of the s-f model. 
Whenever correlation functions that carry a dependence 
on the properties of localized magnetic moments emerged 
we had to perform some crude approximations. It is 
connected to this fact, that we only considered a k- 
independent self-energy by taking the average over the 
whole Brillouin zone. In this direction their is certainly 
room for further improvements. 
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APPENDIX A: HIGH-ENERGY EXPANSION 

The high-energy expansion of the Green's function can 
be obtained from its representation as a resolvent Q): 



mi' 



'k(T 



1=0 



1=0 



l+l 



(Al) 



The coefficients mfl are called spectral moments and are 
determined by an Z-fold commutator with the Hamilto- 
nian Ti. Apart from the mean-field contribution to the 
frequency matrix (0), the self-energy is identical to the 
memory matrix ((7J). Hence, its high-energy expansion 

Sk. = E C^k"Vw" is given by 



{QCcl^\iQCQr\QCcl^ 



,rn+l 



(A2) 

Due to the properties of Q, the coefficients C^™'' can be 
expressed in terms of spectral moments: 



C, 



(1) (2) 



.(1) 



(2) (3) o (2) (1) 



*ko- 



(A3) 
(A4) 



However, in our second-order perturbation theory the Li- 
ouville operator C is replaced by a part C^: 



^(SOPT) 
■"kcr 



(A5) 



k(T 



QCcl 



k(T 



,m+l 



m=0 



If this result is used for the MPT the high -energy 
expansion of the self-energy is given by: 



^(MPT) 

^k<T 



(A6) 

A comparison with the exact expression (|A2I) shows, that 
7k(T is correct to order cj"^. To ensure correctness to 
the same order for "^^^"^^ the parameter akcr has to be 

chosen as 1. To order the self-energy sj^'^^"^' has the 
coefficient 



QCcl 



kCT 



(A7) 



In order to ensure that also this coefficient is exact, it 
has to be equal to 



(A8) 



Equality can be obtained if 6ko- is chosen as suggested 
in equation (|18|l above. As can be seen from (|A4|) in 
the high-energy this order expansion implies the correct- 
ness of the four moments , ■ • • , m!^^ of the Green's 
function. 



APPENDIX B: ABBREVIATIONS 

For the sake of brevity we have introduced some short- 
hand notations in this papers. The full expressions are 
given here. 

Equation ifTOjl : 

= z„{S''n^„) -{S-''clc_^) (Bl) 
qa - 42^(S'^)p^-2(a;^-p^) (B2) 
x„ = S{S +l){n^) - z„{S'n^) + 2z^{S''n„h^^) 

+ {{S^' f{n_^ - K)) - Is'^S^cl^c. + h.c.) 

+ h.c.\ (B3) 



We have evaluated Xa by making use of the equivalence 

= -J(a;, -5(5+1) (n„,)-p,) (B4) 

and arguing that the left hand side vanishes for almost all 
parameter configurations and in particular if mean-field 
decoupling is applied. 
Equation ifT^ : 



X = JZ- 



2J 



s{s + i)-z„{s^) - {s^y 

~z^ {S^){{S^f)-z. (S^f 



S{S+l)-z,{S-)-{S^y E-To-Jz,{S-) 



2JZz„(S^) 



((5^)2) - {sr 



Y = Z- 



s{s + i)-z,{s^) - {s^y 



1 - J 



E-To-.Jz^{S') 



X 



z = 



1 



E -Tq + Jz„{S^) E -To + Jz-aiS"- 



(B5) 
(B6) 
(B7) 



(5^) vanishes in the paramagnetic regime. Hence, Z be- 
comes zero, which implies the same for X = Y = Q. 
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